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Abstract. The thermal evolution of the intergalactic medium (IGM) 
depends on the reionization history of the universe. Numerical simula- 
tions indicate that the low density IGM, which is responsible for the low 
column density Lya forest, follows a well defined temperature-density re- 
lation. This results in a cut-off in the distribution of line widths as a 
function of column density. We use hydrodynamic simulations to cali- 
brate the relation between the cut-off and the temperature-density rela- 
tion and apply this relation to Keck spectra spanning a redshift range 
z = 2-4.5. We find that the temperature peaks at z ~ 3 and interpret 
this as evidence for reheating due to the reionization of helium. 



1. Introduction 

Quasars have provided us with a unique probe of the high redshift universe. 
These bright point sources shine like a flashlight through space, revealing the 
presence of baryonic matter through the light it absorbs. Thus every quasar 
spectrum contains a one-dimensional map of the distribution of matter along 
the line of sight. The extraordinary quality of the spectra obtainable with the 
HIRES spectrograph on the Keck telescope, enables us to extract the wealth 
of information that has been collected by the quasar's light along its journey 
through space and time. Computer simulations of structure formation have 
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been remarkably successful in reproducing these observations. They show that 
the physics governing the high redshift intergalactic medium (IGM), which is 
responsible for the low column density absorption lines (the so-called Lya for- 
est) is relatively simple. The IGM, which contains most of the baryons in the 
universe, is photoionized and photoheated by the collective UV radiation from 
young stars and quasars. On large scales its dynamics are determined by the 
gravitational field of the dark matter, while on small scales gas pressure is im- 
portant. The availability of superb data and a detailed physical model, have 
made the Lya forest into a powerful probe of the high-redshift universe. 

Since shock heating is unimportant in the low-density IGM, most of the gas 
follows a simple temperature-density relation which is the result of the interplay 
of photoionization heating and adiabatic cooling due to the expansion of the 
universe. For densities around the cosmic mean, this relation is well-described 
by a power-law, T = T${p/ p) 1 ^ 1 (Hui & Gnedin 1997). At reionization the gas 
is reheated, resulting in an increase in To and a decrease in 7 (provided that the 
gas is reionized on a timescale short compared to the Hubble time) . In ionization 
equilibrium, To decreases and the slope of the effective equation of state steepens 
(i.e. 7 increases). However, because the timescale for recombination is long, the 
gas retains some memory of how and when it was reionized (Miralda-Escude & 
Rees 1994). 

The distribution of line widths (6-parameters) depends on various mecha- 
nisms. Thermal motions of the hydrogen atoms broaden the HI absorption lines 
and other processes, such as the differential Hubble flow across the absorbing 
structure and bulk flows, also contribute to the line widths. However, the mini- 
mum line width is set by the temperature of the gas, which in turn depends on 
the density. A standard way of analyzing Lya forest spectra is by decomposing 
them into a set of Voigt profiles. Since the minimum line width (6-parameter) 
depends on the temperature, and since column density (N) correlates strongly 
with physical density, there is a cut-off in the b(N) distribution which traces 
the effective equation of state of the IGM (Schaye et al. 1999; Ricotti, Gnedin 
& Shull 2000; Bryan & Machacek 2000). We have used this relation to measure 
the thermal evolution of the IGM from a set of nine Lya quasar absorption line 
spectra. 

This work is more fully described and discussed in a forthcoming publication 
(Schaye et al. 2000). 

2. Method 

We have measured the b(N) cut-off for a set of nine high-quality Lya forest 
spectra, spanning the redshift range 2.0-4.5, eight of which were taken with the 
HIRES spectrograph of the Keck telescope. We used hydrodynamic simulations 
to calibrate the relations between the b(N) cut-off and the temperature-density 
relation. Except for the two lowest redshift quasars, the Lya forest spectra were 
split in two in order to take into account the significant redshift evolution (Az ~ 
0.5) and signal-to-noise variation across a single spectrum. The calibration was 
done separately for each half of each observed spectrum. The synthetic spectra 
were processed to give them identical characteristics (resolution, pixel size, noise 
properties, mean absorption) as the real data. The same Voigt profile fitting 
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Figure 1. The evolution of the temperature at the mean density (left 
panel) and the slope of the effective equation of state (right panel). Hor- 
izontal error bars indicate redshift interval spanned by the absorption 
lines, vertical error bars are 1 a errors. Different symbols correspond 
to different quasars. See text for a description of the models. 



package (an automated version of VPFIT (Webb 1987)) was used for both the 
simulated and the observed spectra. 



3. Results and discussion 

The measured evolution of the temperature at the mean density and the slope 
of the effective equation of state are plotted in Figure 1. From z ~ 3.5 to 
z ~ 3.0, Tq increases and the gas becomes close to isothermal (7 ~ 1.0). This 
behavior differs drastically from that predicted by models in which helium is 
fully reionized at higher redshift. For example, the solid lines correspond to 
a simulation that uses a uniform metagalactic UV-background from quasars as 
computed by Haardt & Madau (1996) and which assumes the gas to be optically 
thin. In this simulation, both hydrogen and helium are fully reionized by z ~ 4.5 
and the temperature of the IGM declines slowly as the universe expands. Such a 
model can clearly not account for the peak in the temperature at z ~ 3 (reduced 
X 2 for the solid curves are 6.9 for To and 3.6 for 7). Instead, we associate the 
peak in To and the low value of 7 with reheating due to the second reionization 
of helium (Hell — > He III). This interpretation is supported by measurements of 
the SilV/ClV ratio (Songaila 1998, but see also Boksenberg, Rauch, & Sargent 
1998 and Giroux & Shull 1997) and direct measurements of the Hell opacity 
(Heap et al. 2000 and references therein). 

The dashed lines in Figure 1 are for a model that was designed to fit the 
data (reduced x 2 is 0-24 for To and 1.38 for 7). In this simulation, which has a 
much softer UV-background at high redshift, Hell reionizes at z ~ 3.2. Before 
reionization, when the gas is optically thick to ionizing photons, the mean en- 
ergy per photoionization is much higher than in the optically thin limit (Abel 
& Haehnelt 1999). We have approximated this effect in this simulation by en- 
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hancing the photoheating rates during reionization, so raising the temperature 
of the IGM. 

Since the simulation assumes a uniform ionizing background, the tempera- 
ture has to increase abruptly (i.e. much faster than the gas can recombine) in 
order to make 7 as small as observed. In reality, the low-density gas may be 
reionized by harder photons, which will be the first ionizing photons to escape 
from the dense regions surrounding the sources. This would lead to a larger 
temperature increase in the more dilute, cooler regions, resulting in a decrease 
of 7 even for a more gradual reionization. Furthermore, although reionization 
may proceed fast locally (as in our small simulation box), it may be patchy 
and take some time to complete. Hence the steep temperature jump indicated 
by the dashed line, although compatible with the data, should be regarded as 
illustrative only. The globally averaged To could well increase more gradually 
which would also be consistent with the data. More data at z ^ 3 is needed to 
determine whether the temperature rise is sharp or gradual. On the theoretical 
side, more realistic models should include radiative transfer effects, which are 
important during reionization. 

Together with measurements of the He II opacity, which probe the ionization 
state in the voids, the thermal history of the IGM provides important constraints 
on models of helium reionization. Furthermore, the temperature of the IGM 
before the onset of helium reionization can be used to constrain the redshift of 
hydrogen reionization, which marks the end of the dark ages of cosmic history. 
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